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We argue that the holographic description of four-dimensional BPS black holes naturally 
includes multi-center solutions. This suggests that the holographic dual to the gauge theory 
is not a single AdS2 x S 2 but a coherent ensemble of them. We verify this in a partic- 
ular class of examples, where the two-dimensional Yang-Mills theory gives a holographic 
description of the black holes obtained by branes wrapping Calabi-Yau cycles. Using the 
free fermionic formulation, we show that 0(e~ N ) non-perturbative effects entangle the two 
Fermi surfaces. In an Euclidean description, the wave-function of the multi-center black 
holes gets mapped to the Hartle-Hawking wave-function of baby universes. This provides a 
concrete realization, within string theory, of effects that can be interpreted as the creation 
of baby universes. We find that, at least in the case we study, the baby universes do not 
lead to a loss of quantum coherence, in accord with general arguments. 
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1. Introduction 

The study of quantum aspects of black holes has led to important progress in a deeper 
understanding of quantum gravity. One basic notion is that of black hole entropy, which 
was predicted by Bekenstein |lj and Hawking ||, through semi-classical reasoning, to be 
one quarter of the area of the horizon in Planck units. More recently it was shown in the 
context of string theory || that for special classes of black holes the microstates of the 
black hole consist of bound states of suitable configuration of branes (see JfHU for reviews 
of this subject). In particular it was found that the semi-classical reasoning of Hawking 
agrees with the leading large charge entropy of black hole microstates constructed within 
string theory. 

However, in string theory one can go further and compute, in addition, the subleading 
corrections to the black hole entropy. More specifically, in the context of certain extremal 
black holes obtained in compact ificat ions of type II strings on Calabi-Yau three-folds these 
corrections are captured by topological string amplitudes 0. These results have recently 
led to a concrete formulation of the quantum corrected black hole entropy to all orders 
in string perturbation theory ||. This states that the partition function of a statistical 
ensemble of black hole states Zbh is given by the norm-squared of the topological string 
wave-functionB on the corresponding Calabi-Yau three-fold: 

Zbh = \^top\ 2 - (1-1) 

Since in this relation the string coupling constant is inversely proportional to the charge, 
large charge black holes get mapped to the weak coupling limit of topological strings. 

The appearance of the notion of a wave- function in the above formula, at first sight, 
sounds surprising since we are dealing with a partition function. This was explained in 
||12|| , where it was identified as a Hartle-Hawking wave-function associated to a radial 
quantization, as contrasted to the usual temporal quantization, of the Euclidean black 



hole geometry. In particular, as was observed in ||12|| , the BPS mini-superspace Hilbert 
space TLm (for type IIB compactifications) corresponds to the geometric quantization of the 
phase space H 3 (M), where M is the corresponding Calabi-Yau three-manifold. Moreover, 



as 



1 In this paper, we refer to the topological string partition function ijjto V = exp (^2 g F, 
a wave-function following the interpretation [P,|T0[] of the holomorphic anomaly equations j ffi| for 
V>top as representing the background independence of the geometric quantization of the tangent 
space to the moduli space of the Calabi-Yau three-fold. 
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fixing the electric/magnetic fluxes (Q,P) leads to a distinguished state in this Hilbert 
space 

\Q,P) = e^ XI - plF ^\o,o)en M , 

where X T ,Fi are suitable (canonically conjugate) operators, with the property that the 
black hole state degeneracy Q(Q,P) corresponds to 

n(Q,p) = (Q,p\Q,p}= I n#V«,p)MI 2 - 

J i 

In the real polarization on H (Af), this wave-function is related to the topological string 
partition function via 

Some aspects of the background independence of the black hole entropy and its relation 
to this Hilbert space has been discussed in []T3[ . 



The prediction (|1.1|) has been verified in a number of examples [l|] , 15| , |16 , 17[18 ,lj 
In particular, in the context of a T 2 embedded in the Calabi-Yau, it was shown in ||14| 
that the bound state of D4, D2 and DO branes maps to the partition function of U(N) 
2d Yang-Mills on T 2 . Here the number of D4 branes corresponds to the rank N of the 
gauge group and the chemical potentials for D2 and DO branes can be identified with some 
combination of the theta angle and the gauge coupling of the Yang-Mills theory. The fact 
that in this case the Yang-Mills partition function takes the form of the norm-squared of 



a holomorphic object follows from the results in |20| where the 't Hooft large N limit of 
Yang-Mills theory was studied. 



However, as was noted in ||14j| , there are additional non-perturbative corrections (be- 
having like e~ N ) to the large iV limit which destroy the holomorphic factorization property 
( |1.1| ). The lack of exact factorization is best understood in the free non-relativistic fermion 
formulation of Yang-Mills theory where the two Fermi surfaces are entangled at finite N. 
Our main goal in this paper is to study these corrections and interpret their physical 
meaning in the dual superstring theory!. 



2 These non-perturbative effects also explain the origin of apparent discrepancies in some of 
the examples studied in [21], where computations are done in the strong coupling regime, g s 3> 



1. Because of the non-perturbative corrections, we expect O(l) corrections to the perturbative 



formula (1.1) in this regime and thus we find the question is reversed: Why did some of the 



examples in [ ^T| work at all? 
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What we find is that the correction terms to the large iV limit of the D-brane gauge 



theory can be interpreted as arising from multi-center black holes ||22|| . a special case of 
which correspond to the Brill instantons |[23[ . This, in particular, leads to the statement 
that, while the perturbative 1/N expansion holographically describes a single black hole, 
via its non-perturbative 0(e~ N ) effects, the gauge theory is actually dual to a coherent 
ensemble of black holes. This is an interesting twist to the notion of holography and may 



lead to a resolution to the puzzles raised in in the context of A0IS2 holography. From the 



viewpoint of radial quantization ||12|| , our conclusion is therefore that the Hartle-Hawking 
wave wave-function is not a single universe wave function, but rather the wave-function 
for an ensemble of baby universes. 

In particular, in this language we find that the suitable wave function belongs to the 
"third quantized" Hilbert space, where we have one Hilbert space per baby universe!. We 
find that in order to capture non-perturbative corrections we need to consider the total 
Hilbert space 

n>0 

We find that, if we fix the total flux corresponding to electric/magnetic charges (Q,P), 
the complete quantum state \i/Jtotai(Q, P)) receives contribution from the Hartle-Hawking 
wave- function of an arbitrary number of baby universes and belongs to (Bn'Hfp- It takes 
a particularly simply form. The n universe state is essentially given by 

\MQ,P)) = \Qi,Pi;Q2,P2-,---;Q n ,Pn) 

Qi+...+Q„=Q, 
Pt+...+P n =P 

with a suitable range of sum over charges (and modulo some positivity constraint) and 
where 

|Qi, Pi; Q 2 , P 2 ; ■ ■ ■ ; Qn, P n ) = \Qi, Pi) ® |Q 2 , P 2 ) ® • • • ® \Q n , P n }- 

This multi- universe wave- function is related to the net black hole entropy 0(Q, P) by the 
formula 

00 

Y,(-i) n ~ 1 Cn-i(MQ,P)\MQ,P)) = n(Q,P), 



String theoretic realization of baby universes has recently been discussed in a different context 

in 
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where C n , the n th Catalan number, counts the number of planar binary trees with (n + 1) 
branches, i.e. the number of distinct ways the baby universes can be produced. The sign 
factor (— l) n_1 is perhaps unexpected. This extends the analysis of (T^j to the wave- 
function of universes which are spatially disconnected. 

There is a natural interpretation of this simple factorization structure from the view- 
point of the dual gravity solutions, which turn out to be multi-center black hole solutions 
. Each of the component wave- functions is associated to the near horizon geometry of 



the corresponding black hole. The structure we have found for the multi-baby universe 
wave-function suggests that there is no loss of quantum coherence, in line with the pre- 
dictions of Coleman [P6[ . Essentially, this is because, if we measure the fluxes through 
one baby universe (Qi,Pi), the corresponding wave-function is fixed to be \Qi, Pi) and is 
independent of the other fluxes (Q r , P r ) r ^i or the degrees of freedom on the other universes. 

The organization of this paper is as follows: In section 2 we review the 2d Yang-Mills 
partition function and its large N limit. We also review its relation to black hole entropy. 
In section 3 we consider 0(e~ N ) corrections to the large iV limit of Yang-Mills theory and 
their relation to BPS partition functions of D-brane systems. We give an interpretation 
of these effects in terms of multiple wave-functions. In section 4 we review the relevant 
multi-center gravity solutions. Finally, in section 5 we interpret the non-perturbative large 
N corrections in the context of the holographically dual gravity solutions. 



2. Two-dimensional Yang-Mills and topological string theory 

2.1. Topological strings on local T 2 



Let us briefly review the set-up of Jl4] that relates topological strings, two-dimensional 
Yang-Mills theory and black hole degeneracies. The starting point is a non-compact Calabi- 
Yau manifold M given by the total space of the following rank two vector bundle over a 
two-torus 

0(m) © O(-m) — > T 2 . (2.1) 

Here m is a given integer. This non-compact space can be considered as the local neigh- 
borhood of an elliptic curve embedded in a compact Calabi-Yau manifold. 

We will consider the A-model topological string on the geometry (|2.1] ). The string 
partition function ifjto P (t, g s ) will depend on the cohomology class t G H X,1 {T 2 ) of the 



complexified Kahler form k on T 2 and the string coupling constant g s . It has a perturbative 
expansion of the form 

ipto P {t,g s ) = exp^2g 2 s 9 ~ 2 F g (t), 

g>0 

where F g (t) is the contribution at genus g in string perturbation theory. The stringy con- 
tributions to these perturbative terms can be viewed as generated by world-sheet instanton 
effects. If the Gromov-Witten invariant Nd, g denotes the "number" of instantons of degree 
d and genus g (it is in general a rational number), then 



dt 
d>0 



These contributions are only non-zero for g > 1. In addition to these world-sheet instantons 
effects, there are classical contributions at genus zero and one, given by certain intersection 
numbers. In general for a non-compact target space these are a bit ambiguous, but in the 



case of M these can be computed to be [14 



6 m z 24 

After including these classical parts, the net partition functions F g (t) are quasi-modular 
forms of weight (6<7 — 6) under the usual action of SL(2, Z) on r = it/2n. For example, 



27,28 



F 1 (t) = - log 77, 

F 2 (t) = — - — (lOEi - 6E 2 E 4 - 4£ 6 ) , 
V ; 103680 V 2 ; ' 

with rj the Dedekind eta- function and E n the Eisenstein series, being (quasi) modular 

functions of weight n. This modularity can be understood from applying mirror symmetry 

to the T 2 that turns r into the complex modulus of the dual elliptic curve |?9| and has 



been rigorously proven in [ 30f| - Furthermore, if a suitable anti-holomorphic dependence is 



added, in accordance with the holomorphic anomaly of the string partition function |31J] 
full modularity is restored |52| . 



2.2. D-Branes and 2d Yang-Mills theory 

We will now consider a type IIA compactification on the Calabi-Yau space M. We can 
wrap D4-branes, which we take to cover the base T 2 and one of the two complex fiber 
directions. (This breaks the symmetry m — > — m.) This will give us a 4 + 1 dimensional 
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supersymmetric U(N) gauge theory, that describes a point-particle in the four non-compact 
dimensions. We can further consider bound states with N2 D2-branes, that wrap the T 2 , 
and iVo DO-branes. These lower-dimensional branes will be represented by the Chern 
classes c\(E) and 02(E) that capture the topology of the gauge bundle £ of the D4-brane. 

After taking into account the back-reaction of the supergravity, this collection of D- 
branes will manifest itself as a charged four-dimensional black hole. Since we do not have 
any D6-brane charge, the electric and magnetic charges are respectively given by 

Q = (N 2 ,N ), P=(N,0). 

The black hole partition function that counts the number of BPS states can be identified 
with an index of the corresponding gauge theory. In fact, as explained in || and as we will 
review at greater length in section 4, within the string theory context it is more natural to 
compute this partition function in a mixed ensemble, where we fix the magnetic charges P 
and introduce chemical potentials for the electric charges Q. In this case this means that 
we fix the rank N of the gauge group and sum over the different topologies of the gauge 

bundle. The chemical potentials for DO and D2 branes can be identified with and 

^ 9s g s 

respectively, where the angle 6 is the coefficient of the Tr F A k term in the 4d action. 
So the gauge theory/black hole partition function takes the form 

Z N (g„e)= £ n(W;W„W.).e«p(-^-^). 

N 2 ,N Q >0 \ 9s 9s J 

Here Q(N; N2, Nq) denotes the index of BPS bound states with the given charges. In 
order to relate this black hole partition function with the closed topological string partition 
function as in one needs to make the following identification of the closed and open 

string moduli [|14| 

t = ^g s mN + id. 



As is explained in ||14|| , this particular D4-brane system can be further simplified and 
in fact be identified with a two-dimensional bosonic Yang-Mills theory on the compact 
base T 2 . The action of this model is best written in Hamiltonian form with F the field 
strength of the gauge field A, together with an additional adjoint scalar $ (that is the dual 
momentum to A) 

S = — I Tr ( <$>F + 6$ + -m $ 2 
9s Jt* V 2 
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-N/2 

(a) 



Fig. 1: The spectrum of two-dimensional Yang-Mills is described by a system of 
N free fermions. Depicted are the ground state (a) and small excitations above 
the ground state (b), given by fluctuations of the two Fermi surfaces at p = ±JV/2. 



Here we have set the area of the two-torus to one. Then we can identify the gauge coupling 
with 

9ym =m- g s . 

We can also absorb the factor m by redefining g s — ► mg s , at least if m is non-zeroi. 

Since the solution of two-dimensional Yang-Mills theory on a general Riemann surface 



is well-known [^ , |^J^ , p6| , |37|1 , these identifications give the following exact expression for 
the black hole partition function expressed as a sum over all irreducible representations R 
of U(N) 

Z N = ^ e -|9aC 2 (H)+ieC 1 (K)_ 
R 

Here C±(R) and C^-R) are the first and second Casimirs of the representation R. 

2. 3. Free fermion system 

Two-dimensional Yang-Mills theory on a torus has an elegant reformulation in terms 
of a system of iV non-relativistic free fermions moving on a circle P^JHSd - With natural 
anti-periodic boundary conditions, these fermions have half-integer quantized momenta 



4 Note that this chooses the sign of the string couping g s . Changing m — > — m, which corre- 
sponds to picking the other line bundle to wrap the D4-brane, can be compensated by g s — > — g s - 



In the fermion correspondence, a YM state labeled by a definite irreducible representation 
R of U(N) is given by filling some particular levels pi, . . . ,pjy- The ground state, given 
by the trivial representation, is obtained by filling the states from p = — ^N to p = +^iV, 
see fig. 1(a). The non-trivial representations correspond to the excitations of the top and 
the bottom Fermi levels, as depicted in fig. 1(6). 

The Casimirs of the Yang-Mills representations can be expressed as the total energy 
and momentum of this iV-fermion state 

l -C 2 {R) = E - E Q , C 1 (R)=P, 

with 

N „ N 



i 

Here Eq is the ground state energy 



T 

i=l i=l 



£o = ^(iV 3 -iV). 



It is convenient to add this overall energy shift to the Yang-Mills theory and define 



f ermions 



e -g s E+i8P^ 



This shift \Ci — ► \C 2 + Eq is also natural from the string perspective, since it produces 
exactly the classical contributions to the topological partition function WM 



E o = -FtL(t) - F t (t) 



topK^J top\ Ji 



where 



gj 6gj 24 



with t = ±g s N + id. 

2.4- The large N limit and chiral fermions 
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Fig. 2: In the large N limit the two Fermi surfaces at p = N+ and p = —N- 
decouple and we obtain two independent relativistic free fermion field theories. 



As discovered by Gross and Taylor [2D] , two-dimensional Yang-Mills simplifies consid- 
erably in the large iV limit, defined as 

N — > oo, g s — > 0, with g s N fixed, 

The simplification is most easily understood in the fermionic reformulation. Here the dy- 
namics of the large iV limit is in good approximation described by independent fluctuations 
of the two Fermi surfaces that will be separated by a distance N, see fig. 2. In fact, in 
the U(N) gauge theory (as contrasted with the SU(N) theory) the Fermi levels can be in 
general position at p = N + and p = —N- (upto a shift by \ that we will ignore) as long 
as 

iV + + iV_ = N. 

The partition function now simply factorizes, at least to all order in 1/N asi 

Z N (9,g s )= J2 ^N + (e,g s )-^ N je,g s ). (2.2) 

N + +N_=N 

Here the chiral contribution ifiisr, is captured by the zero fermion number sector of a 
two-dimensional chiral fermionic field theory. It is described by removing and adding 



5 There are non-perturbative contributions that preserve the factorized form, and they have 
been studied recently in |]40|j4l|1 . In the next subsection, we will find other non-perturbative effects 
which break the factorized structure. 
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arbitrary numbers of fermions close to the top Fermi surface. These states have momentum 
p = N + + k with \k\ <C N. Therefore their contribution to the partition function is given 
by 

9sE = ^g s p 2 = const + N + g s k + ^g s k 2 . 

Since the total numbers of particles and holes are equal, the constant contributions cancel. 
And in the large N limit the quadratic term in k is order 1/N and can be ignored in 
the leading approximation. Therefore this sector can be described by a set of relativistic 
fermions with a linear dispersion relation. Including the #-term it is given by tk, with the 
complexified 't Hooft coupling 

t = g s N + + i6. 

The chiral partition function therefore becomes a holomorphic function of t 

ipN+{0,g s ) = il>(t,g a ). 

There is a similar term if; N _(9,g s ) coming from the negative Fermi surface at p = — iV_. 
Here the moduli combine to give an anti- holomorphic function ij)(t, g s ) of 

t = g s N_ - id 

The quadratic corrections to the energy are controlled by g s ~ 1/N and capture the 
perturbative string loop corrections. 

3. Non-perturbative corrections to chiral factorization 

The counting of states of 2d Yang-Mills in terms of the chiral fermions is clearly only 
an approximation, because the two Fermi surfaces are not independent. Particles can 
move from the top to the bottom levels. These order e~ N effects give a non-perturbative 
entanglement of the two chiral systems as we will now explain. 
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3.1. Overcounting fermion configurations 



These corrections can be elegantly described as follows. In the large N limit we are 
treating each of the two Fermi levels as the surface of an infinite deep sea, while in fact 
the sea is only of finite depth N . Therefore we make mistakes if we create holes which lie 
too deep under the surface. 

To be as concrete as possible, let us first decompose the partition function Zn in 
super-selection sectors 

Zn = ^ Z N+;N _ , 

N + +N^=N 

where Z^ + ,N- receives contributions from configurations where there N + fermions with 
positive momenta Pi > and iV_ fermions with negative momenta pi < 0. We can think 
of such a configuration as being created out of a Fermi sea, where the positive level is 
at p = N + and the negative level at p = — iV_, as we did in the previous section. Now 
by convention we will assume that all the holes that are created with positive momentum 
have gone to excitations of the top Fermi surface, and that all the negative momentum 
holes have gone to the bottom Fermi surface excitation. 

If we denote the chiral amplitudes corresponding to this configuration as and 
V'tv. j then naively we can make the approximation 

Zn+,N_ ~ i>N + ^ N _- 

However, this expression is clearly incomplete, since it will overcount states. Each chiral 
wave- function does not limit the momenta of the holes to be positive for ipN+ or negative 
for ip N _ respectively. In ipN+ there are configurations where, for instance, n_ holes have 
been made with negative momenta in order to create n_ particles with positive momentum, 
raising the top level of Fermi sea from N + to N + + n_. Similarly, ip N _ takes into account 
states where n+ positive momentum holes are created that go to states with negative 
momentum, lowering the bottom level to — iV_ — n+). Typical examples of these "wrong" 
configurations are illustrated in fig. 3(a). 
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(a) (b) 

Fig. 3: The chiral product ?Pn + iPn_ overcounts states. As illustrated in (a), it 
contains configurations where n_ negative-momentum holes are brought to the 
top Fermi level and n+ positive momentum holes are brought to the bottom level. 
These states, that look like (b) in the Yang-Mills theory, have already been counted, 
and have therefore to be subtracted from the chiral product. 

In the Yang-Mills theory these configurations correspond to states with a total of 

= iV_|_ + n_ — n_|_ positive momentum fermions and a total of N'_ = iV_ + — n_ 
negative momentum fermions, see fig. 3(6). These are perfectly fine states, but they have 
already been counted in Z^ 1 n' , under the assumption that all positive momentum holes 
went to the top level and vice versa. Therefore these states have to be subtracted from 
the product i^n+^n^- 

Of course, fig. 3(6) only depicts the corresponding ground state, that now has four 
Fermi surfaces. We note that the topology of the Fermi sea has changed. A "bubble" of 
holes is made with the Fermi sea. In general each of these four Fermi surfaces will have 
fluctuations, as is illustrated in fig. 4(a). In the large N limit, and also when n± ^> 1, we 
can describe these fluctuations as essentially independent. 

It is easy to derive in this way an exact recursion relation for the partition function 
Zn in terms of the mistakes made in the chiral approximation. As we have explained 
above, in that approximation we overcount a total number of n = n+ + ri- holes that have 
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(a) 



(b) 



Fig. 4: The fundamental recursion relation: (a) The configurations that we have 
overcounted create a "bubble" of holes in the middle of the Fermi sea. (b) These 
states can be considered as a tensor product of two copies of the system, where 
the first factor consists of the usual particles, but the second factor is made out of 
holes. 

been created on the "wrong side". Since a hole denotes the absence of a particle, their 
partition function is given by 

Z%° les (e,g s ) = Z n (-9, -g s ). 

(We subsequently do not keep track of the 9 dependence, which is identical for holes and 
particles by parity symmetry 9 — > —9, p — > —p.) On top of these holes we have the usual 
particle states. The total number of these particles is now N + n and they are counted 
by the usual partition function ZN +n (g s ). Since the fermions are free, the total partition 
function of particles and holes is simply the tensor product. This is illustrated in fig. 4(6). 
The partition function of these states comes with an extra minus sign, because the states 
are overcounted in the large N chiral limit and therefore have to be subtracted. We thus 
derive the fundamental relation 

N 

z n(9s) = ^^k(g s )^ N -k(9s) - z N+n(g s )Z n (-g s ). (3.1) 

k=0 n>0 

Now as it stands, this formula cannot be really interpreted rigorously. Even though the 
LHS makes perfect sense for finite N, the chiral wave- functions ipN that appear on the RHS 
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are only defined as an asymptotic expansion in 1/N. Therefore also the correction terms 
must be considered as formal objects. This is clear, since we flipped the sign of g s = g\ M . 
This makes the theory ill-definedi. We therefore have to give another, physically sensible 
interpretation of this result. 

We will first sketch a rough argument and will then try to make more precise sense 
of this in the next section using formal generating functions. First of all, using CPT 
invariance, we can interpret the contributions of the holes as the partition function of a 
gauge theory of negative rank U(—n) or equivalently as a supergroup of pure fermionic 
rank 0\n. That is, we can write 

Z n (-g s ) = Z- n (g s ). 
So the relation ( |3.1|) can be written as 

N 

z N(g s ) = ^2^k(g s )^ N - k (9s) - z N+n(g s )Z- n (g s ). (3.2) 

fc=0 n>0 

Now we want to argue that, using a suitable analytic continuation, we can replace the sum 
over n > with a sum over n < 0. In the next section will give an argument using formal 
generating functions. This last step n — > — n will bring the relation to its final form 

N 

z N(g s ) = Y ^k(g s )^ N -k(9s) - z N-n{g s )Z n (g s ). (3.3) 

k=0 n>0 

3.2. Generating functions 

We will now give a more mathematical derivation of these results using generating 
functions. 

3.2.1 Warm-up: a case with linear dispersion relation. 

Let us first consider a more elementary example to explain the main idea. Consider N 
free fermions with a linear dispersion relation E = p. Here the momenta p are half-integer 
and (to make the system stable) are taken to be positive, so p G Z> + \. This system is 
related to the special case m = of the Calabi-Yau geometry ( p.l|) . It is also exactly the 

6 This instability reminds one of Dyson's argument |4^| about the non-analyticity of the QED 
perturbative expansion: sending e 2 — > — e 2 makes charged configurations unstable. 
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spectrum of gauged matrix quantum mechanics with a quadratic potential i.e., the 
N x N matrix harmonic oscillator with action 

S = \f dt Tr((A$) 2 -$ 2 )- 
The partition function is now defined as 

Z N (t) = J2 e~ tE - 

states 

It is most simply written down in a grand canonical ensemble with chemical potential \i. 
With x = e _M and q = e _t , it is given by the generating function 

Z(x;t) = % Nz N(t) = ]J(l + xq p ). 

N>0 p>0 



In this case we have of course also a simple exact expression for Zn H3 

Z N {t) = q N ^ 2 \{{l-q n )- X . 

n=l 

When considered as a large N string theory, there are only perturbative terms at genus 
zero and one. So, ignoring the non-perturbative effects, the perturbative answer in the 
large iV limit is given by (we suppress the dependence on t in the following) 

oo 

^ N :=Z^=q N2 ^\{{l-q-)-\ 

71=1 

The notion of such a perturbative part only makes sense for large N . For example, to take 
an extreme case, Z Q = 1, whereas according to the above definition tjj would be given by 
the very non-trivial expression Hn>o(^ — q™) -1 - 

It is easy to see that the non-perturbative corrections are order e~ tN = q . One 
simply writes 

oo 

Z N = ^ N .]J(l-q N+n ). 

n=l 

These terms can indeed not be ignored for small N. 

We will now derive a recursion relation for the exact partition function Zjy. Consider 
the Jacobi triple product identity (or boson-fermion correspondence) 



n>0 p>0 NEZ 



15 



With the above notation it can be written as 



Z(x) -Z{x- X ) = ^x N i) N . (3.4) 
xez 

Let us also introduce a notation for the RHS, the generating function of the perturbative 
partition functions 

tp(x) = 2_. x N ipN- 
Nez 

Note that here N is allowed to run also over negative integers, but V'-JV = ipN- We repeat 
the caveat that there is a physical significance to i/jn only for large N. 

With this notation we can now write the fundamental relation (|3.4|) even simpler as 



Z(x) ■ Z(xT x ) = ifj(x). (3.5) 

Here Z is the exact expression, if) the perturbative approximation; the relation is exact, 
though. When expanded back into components of fixed rank, it gives the recursion relation 

ZN+kZk = IpN • 

k>0 

Since Zq = 1, we can write this more suggestively as 

Z N = tp N - ZN+kZk, (3.6) 

fc>0 

where the second term on the right hand side denotes the non-perturbative effects. Since 
tp n ~ e~ tn / 2 , the k th correction term in (|3.6|) is exponentially suppressed with respect to 



tpN by a factor e~ tNk . 

Now, if N is large, the leading approximation is Z^ « ift^. So, we can recursively 
expand the terms on the RHS, at least if k is also large, that is, of the order of N. Then 
we have 

Zn ~ 4>n - i^N+k^k + o(tp 3 ). 

k>0 

Only the terms with fc > 1 have a well-defined large iV expansion and can therefore be 
sensibly interpreted. This expansion of the recursion relation ( ft.6|) continues. At the next 
order we have 

Zn ~ tpN - ^2 ^N+ki^k + ^2 i J N+k 1 +k 2 ^k 1 i ! k 2 + OOj 4 ). 

k>0 fei,fc a >0 
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3.2.2 Two-dimensional Yang-Mills 

Now we consider the relevant case of 2d Yang-Mills theory. Here the dispersion relation 
is quadratic E = \p 2 instead of linear, as in the previous case. The quadratic corrections 
give rise to the non-trivial string expansion. 

Here too, it is most practical to consider the fermion system in a grand canonical 
ensemble with a chemical potential \i for the number of fermions N. Then we can write 
a compact generating function for the partition functions for all U(N) theories at once. 
With the notation 

x = e _At , y = e %e , q = e~ 9 % 
we then have an infinite product representation 

oo 

z(x-,e,g s ) = J2^ Nz N(o,9s)= J] (± + xy p q p2/2 )- (3-7) 

N>0 V= — 00 

Because p runs over both positive and negative values, Z is an even function of 9. 

In a similar way, there is a simple expression for the chiral partition function 
that describes the fluctuation of the top Fermi surface at p = N + 3> 1. 

(9, g.) = I ^-x~ N + J] (l + x yW*' 2 ) (l + x^y-^^) . 
J mx p>0 

One easily see from this product formula that the N + dependence of ipN + is given by 

ip N+ (#, g 8 ) = y W q W-&"+ -MO- ^9sN +1 g s ). 

The prefactor gives the leading energy and charge of the ground state. The second factor 
is the perturbative expansion and is only a function of g s and of t = g s N + %9. So we can 
verify using these explicit formulas that, as claimed, 

tpN+{9,g s ) = *P(t,g s ) 

with 

m 9s) = i II (l + x e" V 2/2 ) (l + x^e-KqS/ 2 ) . (3.8) 

J mx p>0 

For the bottom Fermi surface at p = —N- we have similarly (with t = g s N- — i9) 

^(t, g s ) := (9, g s )=f^U{ 1 + x e" V /2 ) (l + V pV2 ) • (3-9) 

J mx p>0 
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Again, these expressions only makes sense at large N, as an expansion in 1/JV. 

We will now repeat some of the previous steps that we performed in the linear case. 
Consider the formal product 

+00 

Z(x;e,g s )-Z(x- 1 ;-9,-g s )= J] (l + x/// 2 ) (l + x" W^ 2 ) • (3-10) 

p— — 00 

Here the second factor can be considered as a partition function for holes. By CPT it is 
obtained by reversing the signs of all potentials 

(P,0,9s) -> (-(JL,-8,-g B ). 

This hole factor is clearly problematic when viewed as a power series in q, since the holes 
can have arbitrary negative energy — \p 2 ■ It diverges badly. At this point we can therefore 
only consider it as a formal expansion in powers of g s . 

Now consider the Laurent expansion in x of the above product. This is best done 
by splitting the product over all p in ( p.lOp in a product over p > and a product over 
p < 0. Then, comparing to the expression (|3.8|) for the chiral wave-function, we see that 
(suppressing the ^-dependence in our notation) 

Z(x;g s )Z(x- 1 ;-g 8 )= £ x N ++ N -^ N+ (g s )^ N _(g s ). (3.11) 

N + ,N_>0 

Now we can formally expand the contribution of the holes as 

Z{x~ x \-g s ) = ^ x~ N Z N (-g s ). 

N>0 

This gives an identity of the form 

JV 

^2 Z N+n (g s )Z n (-g s ) = ^2^k(9s)i>N-ki9s)- 

n>0 k=0 



Using Zo = 1 we thus recover relation (|3.5|) 

N 



Z N (g s ) = ^2^k(g s )^ N -k(9s) - ^2 Z N+n (g 3 )Z n (-g s ). (3.12) 



fe=0 n>l 



Now there is a way to make more sense of the holes factor Z(x\ —0, —g s ). Formally, 
using the identity 

(1 + a" 1 ) =a" 1 (l + a), 
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we can write 

oo 

Z(x- 1 ;-e,-g s )= J] (l + z-W p2/2 ) 

p= — oo 

oo 

x a y b q c J] (l + ajy-V /a )~^;fl,^) 

p= — oo 



Here the a, 6, c are some constants, strictly speaking all infinite. But we could assume that 
they are given by, for example, zeta-function regularization. Since the sum runs over both 
positive and negative p these regulated sums are actually zero. 

In fact, we will turn things around and define Z(—g s ) through this procedure. Note 



that equation (|3.11| ) now becomes 



Z( X] 9,g s ) 2 = x N++N -i>N + (g s W N _(9s 

N + ,N^>0 



or, when written in components, reproduces (|3.3j) 



z n(9s) = ^k(9s)ip N -k(9s) ~ Y Z N _ n (g s )Z n (g s ). (3.13) 



fc=0 n>0 



As in the previous example, we can only trust this formula in the large N limit and for 
the terms with k ^> 1 finite. 

We can iteratively solve (|3.13| ) for Zn in a power series expansion in iftkipN-k as 



Z N = J2(-l) n - 1 C n - 1 Yl iP N i + ---iP N ^ N i_--^ N n, (3.14) 

n=l at|h \-N™+Ni-\ \-N"=N 

where C n in the coefficient is the Catalan number, 

C n = (3-15) 
n\(n + 1)1 

This combinatorial factor arises since the generating function C(x) — Y2^=iCn-i% n for 
the Catalan number obeys the quadratic equation, 

C(x) = x + C(x) 2 , (3.16) 

which has the same structure as ( |3.13| ) with the identification x — > —ipip and C —>■ —Z. 
The Calatan number C n -i is known to be equal to the number of binary bracketings of 
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n letters. Or, put differently, it counts the number of ways to create a planar binary tree 
with n branches. An interpretation of this combinatorial coefficient in the context of the 
gravity dual will be discussed in section 5. 

Note that the expansion ( |3.14| ) is reliable only in the regime where the baby universe 
number k is much larger than the corresponding parent K and much bigger than one, 
K 3> k ^> 1. This gives a hierarchical structure to the various terms in ( |3.14| ), which will 
be important for a gravity interpretation of the coefficients C n -±. 

3. 3. A convenient rewriting of the result 

In this section we recast the above result in a convenient form which will be more 
immediately applicable to our holographic gravitational interpretation. 

Let us recall that the number of BPS degeneracies of D4, D2 and DO branes 
0(A^4, N 2 , Nq) can be computed by Fourier transform of the Yang-Mills answer. In partic- 
ular (with iV 4 = N) 

n { N,N 2 ,N a) =f d (L) i (l-) e ,J^ + 2 -^)z N{g ,,0 ) (3.17) 
«/ \ 9 s J \ Q s J \ Q s 9 s J 



Now, let us define 



27r 2 iVn iV 9 7r#\ , /l 



S I •> 



ipN,N 2 ,N o {0, g s ) = exp I 1 — ) ip I ^-Ng s + id, g t 

V 9 s 9 s J \^ 

which we sometimes also denote as 

ipN,N 2 ,N (0,g s ) = (N,N 2 ,N ,9,g s \iP). 

This is the wave-function of the topological string in the corresponding flux sector. For 
brevity of notation, sometimes we use the following notation, 

\N,N 2 ,N ) = (NN 2 ,No\ijj). 
In other words in this notation we would have 

Tl>N,N*,N (e>9s) = (d,g a \N,N 2 ,N Q ), 
which we hope does not cause confusion. 
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Then keeping the leading all order in the 1/N expansion we can write (|3.17| ) as 
n(N, N 2 , N )= Jd (1) d (jj \iPn,n 2 ,n q (9, g s )\ 2 . 
In other words, to all orders in the 1/N expansion 

Q(iV,iV 2 ,iV ) = (N,N 2 ,N \N,N 2 ,N ). 

Note that the extra sum over k in shifting iV up and down between ift and if) is already 
incorporated by the integral over 9/g 8 (which is best seen by analytic continuation), and 
does not need to be included on top of this. 

Now we come to using the formula (|3.13|) to incorporate effects which are order 



exp(— tN) and smaller^. However the corrections involve more than one copy of ip and 
ijj. This suggests that we try to interpret this as a multi-Hilbert space corrections. 

For simplicity let us start with the first application of the recursion formula to ( |3.13|) . 
The recursion relation leads at this order to (where we suppress the extra shifting of N 
between ip and ip as that would be automatically taken care of by the inverse Fourier 
transform discussed above) 

-\il>N-k{8,9s)\ 2 \M0,9e)\ 2 - 

In order to give this expression a Hilbert space interpretation we introduce additional 
variables which are gotten rid of by additional delta function integrals: 

x \if; N - k (e,g s )\ 2 \MO',g's)\ 2 ^P(+N 2 e/g s + N /g s ). 

Next we write each of the delta functions (taking the periodicities of the chemical 
potential into account) as S(X) = J2 m emX ■, which leads to rewriting the above as 



Q(N, N 2 , N Q ) 2 = - J2 \ 9s, 0', g' s ) \iiN-k,N 2 - m ,N o - P (0, 9 



s)\ 2 



fc>0,m,p ' 



7 Note that this also explains why the 4d black hole interpretation of topological strings is not 
a good starting point for degeneracies of 5d black holes where one may take a small number of 
magnetic branes N ~ 1, because these corrections are of order 1. 
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This result can be summarized in the notation we introduced before as 



Q(N, N 2 , JV ) 2 = - Yl ( N -k,N 2 -m,N - P \N-k,N 2 -m,No-p)- (k,m,p\k,m,p). 

k>0,m,p 

The generalization to arbitrary orders in the recursive relation is also clear. We define 
a state in the sum of arbitrary number of copies of the Hilbert space. We define 

n 

iv>»>= E (g)KA«) 

N i N i N i i=l 

Jv 4' lv 2' 

where each term in the sum is restricted by the condition that 

i i i 

and where all N\ > 0. Our result is then summarized as 

Q(iV,iV 2 ,jV ) = J2(-l) n - 1 Cn-l(lpn\lpn) (3.18) 

n 

We will give an interpretation of this result in section 5, after we discuss some aspects of 
the holographically dual gravity solutions. 

4. Gravity interpretation 

In this section we will consider a type IIB compactification on a Calabi-Yau manifold 
M. Mirror symmetry will relate this to an equivalent description within type IIA, which 
we can use to compare with the gauge theory result in the last section. 

Let us start by describing the standard single center black hole solution at the string 
tree level. Consider a static spherically symmetric metric, 

ds* = -^-dt* + ^ £ (dx^ + dsly, (4.1) 

where r = \x\. The metric ds^ Y °f the Calabi-Yau three-fold can also depend on r, so we 
have a "warped compactification." Let X 1 (I = 0, 1, h 2,1 ) be the projective coordinates 
of the complex structure moduli space of the Calabi-Yau manifold normalized as 

S = -^Imr.j X'X J = |lm [X'djFo] , (4.2) 
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where F (X) is the genus topological string partition function and 17 j = djdjF Q is the 
period matrix of the Calabi-Yau three-fold. Assuming that X 1 are also functions of r, the 
classical BPS black hole solution with electric charge Qi and magnetic charge P 1 is given 
by 



Re X 1 = 


P 1 


+ C 1 






\x\ 

Ql 




(4.3) 


Re d T F Q = 


+ d I: 




\x\ 







where c 1 , di are integration constants that correspond to values of the Calabi-Yau moduli 
at spatial infinity. This shows how the Calabi-Yau metric ds^y depends on r = \x\. (The 
Kahler moduli are kept constant.) The horizon is at x — > 0, where X 1 approaches the 
attractor point J0J : 

Re X 7 ~ ^, Re d/F ~ (4.4) 

\x\ \x\ 

independently of their values at the infinity. The semi-classical entropy for this solution is 
given by 

S<g ) H (P,Q) = ^lm[X I d I F ] 

= F {X) + F (X) + tv ^ Qi Im X 1 , (4-5) 

i 

where X Ij s are fixed by ([4.4| ). This is also related to the asymptotic behavior of S(x), 



at \ S gH {P,Q) 

S[x) ~ Dn , x 0, 



and for this reason we may regard S(x) as a "position dependent entropy". The near 
horizon geometry is AdS2 X S 2 x M, where the complex structure moduli of the Calabi- 
Yau three-fold M are fixed by the attractor equations (|4.4| ). 

String loop effects modify the low energy effective action, and the black hole geometry 
is also changed accordingly. In , it was shown how to systematically take these perturba- 
tive effects into account. Remarkably the entropy formula to all order in the perturbative 
expansion can be concisely expressed in terms of the topological string partition function 
Ftop{X), 

oo 

F top {X)=Y,F 9 (X), (4.6) 

3=0 
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where F g is the genus g partition function. In || it was shown that the string loop corrected 
entropy Sbh(P, Q) of [0 can be expressed simply as 



Sbh(P, Q) = F top (X) + F top (X) + Qi&, 



(4.7) 



where the Calabi-Yau moduli are fixed to be 



X 1 = P T + 



(4.8) 



and 



s are non-linearly related to the charges (P, Q) by 



Qi 







(X) + F top (X)] . 



(4.9) 



This generalizes the classical entropy formula ( |4. 5| ) , and a quantum version of the attractor 



equation is given by ( |4.8|) and Q4.9|). One can regard ( f4.7|) as the Legendre transformation 
between the entropy, which depends on P and Q, and the topological string partition func- 
tion Fto P + F top , which naturally is a function of P and Motivated by this observation, 
and the earlier work 0, it was conjectured in that the number 0(P, Q) of microscopic 
BPS states of the black hole is given by the Laplace transformation of the topological 
string partition function as 



and X 7 's are given by (f4~8[). 

The main purpose of this paper is to understand non-perturbative corrections to this 
formula. In this section, we will discuss sources of such corrections from the gravity point 
of view, and compare them with the results we found in the gauge theory point of view in 
the last section. 




(4.10) 



where 



iJto P (X) = expF top (X), 



(4.11) 
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4-1. Multi-Center Solutions 

The crucial observation is that the spherically symmetric geometry given by ( |4.1| )- 
( fl.3| ) is not the only solution for a given set of charges (P, Q) that preserves half of the 
supersymmetry. In fact there are multi-center solutions satisfying the same asymptotic 
behavior at spatial infinity [^2j . Each of these solutions is characterized by a decomposition 
of the charges (P, Q) as 

n n 

P' = ^P/, Q, = X^> (4.12) 

i=i i=i 

where n is the number of disjoint horizons of the solution and (Pi, Qi) are charges associated 
to each horizon (i = 1, . . . , n). 

Before describing general multi-center solutions, it would be instructive to discuss the 
simplest case when there is only one gauge field and no scalar field. (This corresponds 
to the type IIB compactification on a rigid Calabi-Yau manifold, h 2 ' 1 = 0.) Furthermore 
suppose that the black hole carries no magnetic charge, i.e. P = 0. Multi-center solutions 
in this case have been known for a long time, and they are called conformastatic solutions 
|[45|| . When there are several extremal black holes whose charges are of equal sign, their 
gravitational attraction is balanced by the Coulomb repulsion, and they can be placed at 
arbitrary positions in three spatial dimensions and still remain static. Such a solution can 
be constructed as follows. Consider a scalar function S(x) of the form 




S(x)=* c+> — , (4.13) 



whose square- root solves the Laplace equation A^S(x) = in three dimensions. We 
choose Qi to be all positive and the constant c is also positive, so that S(x) never vanishes. 
The metric is then given by 

ds 2 = --^-dt 2 + ^ y (dx a ) 2 + ds 2 CY . (4.14) 

S ^ 71 a=r 2 ,3 



The function a/ S(x) also serves as the scalar potential for the gauge field. This describes 
a collection of extremal black holes with charge Qi at Xi. In fact, as x approaches Xi, S(x) 
behaves as 



S(x) 



12 ' 
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where the geometry approaches that of AdS% x S 2 with charge Qi. Toward spatial infinity, 
we have 

S(x) ~ 7T + -pj- + • • • J , |x| — > oo, 

where 

n 

1=1 

and therefore the solution has the same asymptotic behavior as the single center solution 
with the charge Q. These multi-center solutions also preserve one half of the supersym- 
metry. 



It was pointed out in |M6J that, in the limit of c — > 0, Wick rotation of the solution 



given by ( f4.13| ) and (|4.14[ ) describes quantum tunneling of a single universe into several 



disjoint universes. After the Wick rotation 

TV 

S(x) 77 



ds\ = ^-dr 2 + M (dxr + ds 2 CY , (4.15) 



a=l,2,3 



it is convenient to think of S(x) as the Euclidean time. For S — > 0, the geometry approaches 
to AdS2 x S 2 with charge Q. On the other hand, as S — > oo, the geometry fragmented into 
pieces located at Xi {% = 1, . . . , n), each of which has the form of AdSz x S 2 with charges 
Qi. The Euclidean action / of this solution was evaluated in Hfl as 



I = 7rJ2QiQr ( 4 - 16 ) 



i<j 

It was noted in there that this can be written as 



I = ~\ (s { bUQi) + ■■■ + sM(Q n ) - S^(Q)) > (4-17) 

where S^lj is the semi-classical entropy given by (fi~5|). This suggests [^6| , [24| that e _/ is 
an instanton amplitude for a tunneling of a single universe with charge Q into n universes 
with charges Qi, ■■■,Q n so that its square gives the transition probability expected from 
the detailed balance argument. We will find that the comparison with the gauge theory 
computation in the previous section supports this interpretation. 

This construction can be generalized to the case with h 2,1 > 1 as follows. There are 
(h 2,1 + 1) gauge fields and h 2,1 scalar fields describing the complex structure of the Calabi- 
Yau manifold. First consider the case when the charge vectors (P,,Qj) are all parallel 



26 



to each other. In this case, we are in practice turning on only one linear combination of 
the gauge fields. Since the attractor fixed point depends only on ratios of charges, the 
complex structure of the Calabi-Yau are kept constant. In this case, the horizons can still 
be located at arbitrary points. 



The situation is more subtle with non-parallel charges |2^,^7| . In this case, Calabi-Yau 
moduli at each horizon can be different, and it costs kinetic energy for the scalar fields to 
interpolate between the different horizon values. Moreover the electromagnetic interaction 
does not completely balance the gravitational attraction. The metric in this case takes the 
form 

ds 2 = -^(dt + u;dx 1 ) 2 + V (dx 1 ) 2 + ds 2 CY . (4.18) 



S(x) 



n 



i=l,2,3 

As before, the function S(x) and the Calabi-Yau moduli are combined into (h 2,1 + 1) 
variables X (x) normalized as in (|4.2|) . Defining Qi(x) and P 1 (x) by 



X 1 are determined by 



p 



+ c J 



1=1 



Re X 1 = P'(x), 
Re d T F = Qi(x). 



(4.19) 



(4.20) 



Thus the Calabi-Yau metric ds^y also depends on x. The function S(x) is determined by 



7T 



S(x) = -Im [X I d I F ] 



The off-diagonal component in the metric is found by solving 



tdcu = ( pI ( x ) dQi(x) -Qi dP\x)) 



(4.21) 



where * is the Hodge star operator with respect to the flat metric on R 3 . The integrability 
of (|4.21| ) leads to a set of constraints on the locations of the horizons: 



. , . \*^i dy 7 



^ Re (Qij c 1 -Pi d T ) 



i = l,...,n, 



(4.22) 
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Fig. 5: The metric ( 4.18 ) is sketched in (a), it describes an asymptotically Bat 
space-time with multiple black holes with charges (P/ ,Qn) (i = l,...,n). In the 
near horizon limit (6) this geometry gives a disjoint sum of n AdS2 X S 2 geometries. 



where 



en = Yl (Qn P j - P i Qji) ■ ( 4 -23) 
I 



It should be noted that, in addition to this constraint, we need to require that the factor 
S(x) in the metric remains non-zero. With ( |4.2| ), (|4.19 ) and ( |4.20| ), requiring this for all 



x E R 3 implies an inequality on the charges (P/, Qn). We will see that, in the example we 
studied in the last section, this constraint agrees with what we found in the gauge theory 
side. 

As in the single charge case discussed at ( f4.15|) , the Wick rotation — dt 2 — > +dr 2 of 
( [4. 18| ) gives a metric which is asymptotically flat at spatial infinity and fragments into 
several AdS2 throats for S — > oo (see fig. 5). Because of the off-diagonal term Uidx x dt, 
the metric in general becomes complex-valued after the Wick rotation. Since it becomes 
real in the asymptotic regions, it is still appropriate to use it as a saddle point and it can 



contribute to the functional integral [48 



4-2. Baby universes interpretation 

How should we interpret the existence of such multi-center solutions? It is reasonable 
to expect that the Euclidean functional integral contains a sum over these solutions in the 
saddle point approximation. In the standard AdS/CFT correspondence, the low energy 
limit on D branes is dual to the near horizon limit in the gravity side. For the multi-center 
metric, a natural near horizon limit gives n disjoint universes, each with the AdS2 x S 2 xM 
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S(x) 




Fig. 6: As a function of Euclidean time, which can be identified with the local 
entropy S(x) of the metric ( 4.16[ ), the geometry describes the branching off of baby 



universes. 



S(x) 




Fig. 7: As S 
hole centers. 



oo, the number of branches increases to n, the number of black 



geometry, and the moduli of the Calabi-Yau three-fold M at each universe are fixed by the 
attractor equations. Thus, we are led to conjecture that the full partition function of the 
gravity theory contains a sum over disjoint universes (baby universes), see fig. 5. 

Note that if we consider the projection of our geometry to the Saline, and consider the 
Euclidean space which is the pre-image of the Euclidean time S, we obtain the topology 
of a branched tree, where £ can be viewed as the 'height' of the tree, as shown fig. 6. The 
number of components of the space for a given value of S denotes the number of branches 
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of the tree at a fixed height. As S — > oo the number of branches equals n, the number of 
black hole centers, see fig. 7. 

In fact with a fixed asymptotic flux, we are forced to do this summation since there 
are instantons which change the number of baby universes [f46||2"4|1 as we saw earlier. Since 
each AdS2 throat contributes to the partition function by a factor of exp Sbh(Qi, Pi) 
(i = 1, n), the total partition function should be of the form, 

oo 

Z BH = J2C n -i exp(S BH (Q 1 ,P 1 ) + --- + S BH (Q n ,P n )). (4.24) 

n=l p=p 1 + -.-+p n 

Q=Qi+-+Q n 

The sum over charges may be restricted by the condition that the gravity solution exists. 
The coefficient C n _i is the Catalan number and reflects the distinct planar trees (where 
on each node of the tree the parent universe splits off say to the left) with n branches 
which lead to this baby universe sum. This would be a valid description at least in the 
regime where we have a hierarchical splitting of parent universes to baby universes, i.e. 
when the charges of the parent universe are much larger than those of the baby universes 
as discussed at the end of section 3.2. 

We expect a holographic dual to the gravity theory to capture this branching process, 
and we should be able to interpret the different contributions in the gauge theory as coming 
from the different geometries that can contribute to the Euclidean functional integral. 
Moreover, since the latter have finite action, this translates into gauge theory configurations 
which are weighted by e~ N with respect to the vacuum configuration. In the next section, 
we will show that the e~ N terms in the gauge theory discussed in the previous section 
indeed correspond to baby universes and have the right hierarchical branching structure of 
Q4.24p . Interestingly, we will find that the gauge theory predicts an additional sign factor 
of (-1)™" 1 in (Of. 



5. Comparison between gauge theory and gravity theory 

In the previous section we have argued that the D-brane gauge theory is holographi- 
cally dual not to a single black hole, but to an arbitrary ensemble of extremal black holes 
with a given total flux. In other words, fixing the flux at infinity allows more than a single 
black hole solution. This was used to suggest that the gauge theory is holographically 
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dual not to a single AdS-z x S 2 but to an ensemble of A0IS2 x S 2, s. For each horizon the 
arguments of || apply and 

expS(P,Q) = {P,Q\P,Q), 

where 

(<&|P, Q) = e-^-*Vtop (P + -$) = ^p )0 ($). 

7T 

Using this we can rewrite the relation ([4.24 ) as 



Q(P,Q) = X^n-l (Pi,Qi\Pi,Qi) (5.1) 
n=l Pi.Qi i=l 

where Pi = P and Qi = Q and with the restriction that the gravity solutions should 



exist. This prediction can also be stated using the picture proposed in [|12| for the Hartle- 
Hawking wave-function in the mini-superspace. The case considered there involved type 
IIB string theory on 9-dimensional space being 

X = S 1 x S 2 x M, 

where M is a Calabi-Yau threefold, and the fluxes (P, Q) go through S 2 and some three 
cycles of Calabi-Yau. It was argued that in this case V'p ! q($) is the Hartle-Hawking wave- 
function in the Hilbert space TCm obtained by quantization of H 3 (M) with respect to its 
symplectic structure. 

We now extend this question as follows: Suppose we consider n disconnected copies of 
X, labeled by Xi, and let the flux (Pj, Qi) pierce through each. Then we ask which state 
do we get in Hfj 1 by doing the path-integral on geometries whose boundary is n copies of 
XI From the geometry ( |4.18| ), it is clear that now we get the resulting state 



l^n) =Q$\Pi,Qi)- 

i=l 

This is because the states are ground states of the theory and are determined by the long- 
time evolution, which is precisely the near- horizon geometry of each throat. However we 
can ask a further refined question: Can this ensemble of baby universes be dual to a single 
gauge theory? If this were the case there should be 10 dimensional solutions which connect 
up all the Xi and moreover bound by the constraint that Yli Pi = Pi T^i Qi — Q where 
the P, Q are fixed by the total flux of the brane where the gauge theory lives. These are 



precisely the solutions of [22] discussed above. So we would be instructed to write the 
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Hartle-Hawking state as a sum over all allowed fluxes consistent with the fixed total flux 
and with the constraint that the gravity solutions exist. Let us now see if this expectation 
agrees with the results for the case we have studied, namely the case of the T 2 embedded 
in the Calabi-Yau. In this case we found ( |3.18| ) 

n 

M>»> = E (g)l^,iv*,iv*>, 

Ni,Ni,N* *=1 
oo 

n=l 

where each term in the sum is restricted by the condition that 

i i i 

and where all N\ > 0. This is exactly the same structure anticipated in ( |5.1|) from the 
holographically dual gravity solutions, modulo the factor (— l) n_1 in the inner product, 
which would be interesting to explain from the gravity side. The constraint that N\ all 
have the same sign is clear from the restriction that the gravity solution exist. Note that, 
if the sign of N± changes, then there would be points in R 3 where the "position dependent 
magnetic charges" P 1 (x), discussed in the last section, all vanish. This is because if we go 
from a throat with a positive value of to one with negative value, we will cross a point 
where the D4 brane number is zero. Since there are no other magnetic charges in this 
case, this leads to zero classical entropy for S(x) for some x, and thus the gravity solution 
would become singular. 

5.1. Loss of quantum coherence? 

It is a natural question to ask if our multi-universe Hartle-Hawking state \ift n ) leads to 
a loss of quantum coherence. In a naive sense one may think that it does, in that we have 
a sum over all n-state wave-functions with the total flux condition satisfied. However, as 
it stands the mixture with the other universes is very simple and captured just by some 
global flux conservation. In particular if we measure the flux P ouri Q OU r in our universe the 
wave- function ibp n is completely determined. This is consistent with the proposal 

/ J o it t 7 our A. J xx 

of Coleman [^SJ (see also P9| , |50| , pl| , |52[] ) that the creation of baby universes does not lead 
to a loss of quantum coherence. More precisely it was argued that once one measure the 
coupling constants in our universe we will have a pure state. This is consistent with our 
scenario where measuring the flux in our universe is sufficient to lead to a pure state. 
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5.2. Lessons for Holography 

However our finding raises a more interesting question in the context of holography: 
It has been argued that the existence of a unitary gauge theory which is holographically 
dual to a black hole must automatically lead to a resolution of information puzzle for black 
holes. Here, in the context of our simple example, we are finding that the gauge theory is 
not dual to a single black hole but to an ensemble of them. In such a scenario the unitarity 
of the gauge theory evolution operator may not have a direct implication for the unitarity 
of the physics in a given black hole sector. 

There is a precursor for such a sum over geometries being reflected in a gauge theory. 
This is the case of the finite temperature Yang-Mills theory, where one expects contribu- 
tions from both the thermal AdS geometry as well as the AdS Schwarzschild geometry. 
In that case there is, moreover, a phase transition in the semi-classical limit as one varies 
parameters. This Hawking-Page phase transition of geometries exchanging dominance 
translates into a large N phase transition in the gauge theory []55| . (See also H53|j55| , |5S| , |57[ 
for recent studies of this system.) 

In the case of the 2d Yang-Mills theory on the torus, while we have seen multiple 
geometries contributing, we do not have a phase transition as a function of the couplings 
(or chemical potentials on the gravity side). This is consistent with the expectation on the 
gravity side as well where the single black hole is always the entropically favored geometry. 
It would be interesting to consider the issue of scattering off of these extremal black holes 
and see how the unitarity of the S-matrix for the gauge theory is reflected on the gravity 
side for configurations which can fluctuate to the multi-centered black holes. 

There are other gravity backgrounds whose holographic dual descriptions involve 
free fermion Fock spaces; description of black holes in terms of fundamental strings 
5^ , |59| , |60|j6l|1 , Mathur's picture of black hole quantum states [^2|, BPS states for type 



IIB string on AdSs x S 5 |6^,0, and two-dimensional string theory [p5| , |55|J5^ 1 , etc. It may 
be possible to extend the description of baby universes we developed in this paper to these 
and other cases and learn more about non-perturbative phenomena in quantum gravity 
and string theory. 
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